clear all; close all; clc;
%������ ��������� �������
R = 120;%��
C =33e-12/2;%��
RC=R*C;%���������� ������� ����
%��� �������
Tmod = 10*RC;
T = RC / 1000;
t = 0:T:Tmod;

% �������� �� �������� �����������.������ ������
E  = 0*(t<Tmod/4) + 2*(t>=Tmod/4);

%�������� ������
U = nan(1, length(t));
dU = nan(1, length(t));

%������ ��������� ��������� 
U(1) = E(1); dU(1) = 0;



for k = 2:length(t)
%��� �� ��������� ���������� ����������� ����������
    U(k) = (RC*U(k-1) + E(k)*T) / (RC + T);
 end
  U=U/2;   
figure;
plot(t, E, t, U)
xlabel('t, sec');
ylabel('U, V');
grid on;